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oo ! Abstract 



We summarise existing experimental and theoretical knowledge on the structure 
function F2 in the region of low and low x. The constraints on the behaviour 
■ of structure functions in the limit of = are listed. Phenomenological low 

parametrisations of the structure functions are collected and their dynamical 
content is discussed. The high energy photoproduction and nuclear shadowing are 



00 

o 



also briefly described. Recent update of the low Q , low x experimental data is 



■ given. 

Oh; 

^ 1 Introduction 

> 

^ ! Interest in the low Q^, low x phenomena in the inelastic lepton-hadron scattering is 

^ I connected with experimental constraints that the low x region in the present fixed target 

experiments can only be reached on the expense of lowering the values down to 1 
GeV^ or less (see e.g. ref. 0). Unified treatment of small- and large regions might 
also be of high practical importance for large data analysis. This stems from the fact 
that in all leptoproduction experiments, including those at HERA, a radiative corrections 
procedure has to be applied in order to extract the structure functions from the data. The 
radiative "tails" originating from processes at values from the interval Q^eas ^ ^ 
contribute to measurements at = Q'^e.as ^^nd the knowledge of the structure functions 
in this interval is necessary for the iterative data unfolding procedure 0, 0. 

Remembering that due to the conservation of the electromagnetic current the structure 
function F2 must vanish in the limit ^ 0, the Bjorken scaling which holds approx- 
imately at high cannot be a valid concept at low Q^. Theoretical models assuring 
a smooth transition from the scaling to non-scaling regions and applicable in a wide 
interval are thus necessary for understanding the experimental data and the underlying 
dynamics. The continuity of the physical processes occuring when passing from low- to 
high is illustrated in fig.l. We shall be predominantly concerned with the region of 
low and high W , i.e. beyond the resonances. 
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The need for modifying the QCD improved parton model by including contributions to 
the structure functions behaving as l/{Q'^)'^{n > 1) is clearly visible in the data [Q. These 
contributions, called "higher twists" are important at moderate values of GeV^). 
They follow from the operator product expansion [^] and describe effects of the struck 
parton's interaction with target remnants thus reflecting confinement effects. Clearly the 
theoretical description unifying the confinement and deep inelastic regions should contain 
terms of well defined physical origin, corresponding to such contributions. 

The purpose of this paper is to collect the existing knowledge about the region of low 
Q^, low X (i.e. below 3 GeV^ and x below 0.03 or so) in a way which could at the same 
time make it easy to use in practical applications. The paper should be considered as an 
extension of the comprehensive review of the small x physics towards the more detailed 
treatment of the low problems. We shall be exclusively concerned with charged lepton 
inelastic scattering. The recent review of problems specific for inelastic neutrino (and 
antineutrino) interactions is presented in the ref. [0. It is predominantly the structure 
function F2 which will be discussed; particular final state structures, like jets, diffractive 
dissociation, etc. will not be considered. 

The content of the paper is as follows. After basic definitions and constraints (section 
2) we present theoretical ideas and models which describe the low physics (section 3). 
High energy photoproduction (section 4) is then followed by a description of phenomeno- 
logical parametrisations of structure functions (section 5). Special attention is given to 
dynamical models of a low behaviour of F2 (section 6). Nuclear shadowing is described 
in section 7 and finally an update of experimental data is given in section 8. Section 9 
contains conclusions and outlook. 



2 Basic definitions and constraints 



Kinematics of inelastic charged lepton scattering is defined in fig.2. One photon exchange 
approximation is assumed throughout this paper. The imaginary part of the forward 
Compton scattering amplitude of the virtual photon is defined by the tensor iy^'^(see 




In this equation = — where is the square of the four-momentum transfer , x = 
Q^/ {2p ■ q) the Bjorken scaling variable and M is taken as the proton mass. The invariant 
quantity p ■ g is related to the energy transfer v in the target rest frame, p ■ q = Mu. The 
invariant mass of the electroproduced hadronic system, W, is then W"^ = M"^ + 2Mv — Q"^ . 
Often one denotes W"^ = s. 

The deep inelastic regime is defined as a region where both and 2Mz/ are large and 
their ratio, x, is kept fixed. At smaller than few GeV^, x can probably no longer be 
interpreted as a momentum of a struck parton but it remains a convenient variable for 
displaying the data. The functions Fi{x,Q^) and F2{x,Q'^) are the structure functions 
of the target. For a nuclear target it will be assumed that the structure functions are 
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normalised to the number of nucleons in the target nucleus and denoted F^, i =1,2 
(except for the deuteron where a symbol will be used). The tensor W^^ satisfies the 
current conservation constraints: 



q^.W^" = 
q^Wf"" = 



(2) 



This follows from the fact that W'^'^ is related to the matrix element of the product of the 
electromagnetic current operators (x) : 



(3) 



where the symbol 'T' denotes time ordering. 

Let us rearrange eq.(|l]) in order to display explicitly the potential kinematical singu- 
larities of the tensor W^" at = 0: 



W''-(p,q) 



' Fi ^ F2P- q\ qf^q" F2 p^'q" + p^'q^' 



M M q^ 



q^ 



M 



q^ 



(4) 



These singularities cannot be real and appear only as artifacts of the way we wrote up 
W'^'^. In order to eliminate them we have to impose the following conditions on the 
structure functions Fi in the limit 0: 



Fi ^F2p-q 



(5) 
(6) 



M ' M g2 

These conditions have to be fulfilled for arbitrary u. They will play important role for 
the parametrisations of the structure functions at low Q^. 

The differential electroproduction cross section is expressed in the following way by 
the structure functions F,-: 



dQ^dx ^ 



Mxy\ F2{x,Q^) 



2E 



X 



+ 1 



2m^ 



y'Fi{x,Q') 



(7) 



where E denotes the energy of the incident lepton in the target rest frame, m is the 
electron (muon) mass, y = u/E and a is the electromagnetic coupling constant. 

Instead of Fi the structure function R{x, Q"^) defined as 



i?(x,Q2) = ^ 



;i + 4MV/g2)F2 

2xFi 



- 1 



is often used where ctl and cr^ denote the cross sections for the longitudinally and transver- 
sally polarised virtual photons respectively. The differential cross section (0) then reads 



dQ'^dx 



X 



Mxy 

1 — V h 

^ 2E 



2rrf \ 1/2(1 + 4MV/g2 
"O^j 2{1 + R) 



(9) 
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Real photons are only transversally polarised and therefore and R vanish when — > 0. 
This vanishing follows from eqs (^) and (j^). The function R is related to the frequently 
used longitudinal structure function Q^) via R = Fl/2xFi, and 

F,{x, Q') = [ 1 + j F2 - 2xF,. (10) 

At large Q^, Fl is directly sensitive to the gluon distribution function, which plays a 
crucial role in the interactions at small x, f^. 



3 Basic theoretical concepts relevant for the small 
region 

In the leading logQ^ approximation of perturbative QCD which is applicable in the high 
region the structure function F2{x, Q^) is directly related to the quark- and antiquark 
momentum distributions, qii^x^Q"^) and qi^XjQ"^): 



F,{x, Q')=xY. [g,(x, Q^) + qi{x, Q^)] (11) 



denotes quark flavours and are the quark charges. At high the quark and 
antiquark distributions exhibit the approximate Bjorken scaling mildly violated by the 
QCD logarithmic corrections. The evolution of these distributions with is described 
by the Altarelli-Parisi equations |p. These equations as well as the relation (|Tl]) acquire 
corrections proportional to a;s(Q^) in the next-to- leading logQ^ approximation. 

Systematic analysis of the structure functions in the Bjorken limit can be done using 
the operator product expansion of the electromagnetic currents (cf. eq. ^,0. This 
expansion leads to the expansion of the structure functions in the inverse powers of Q^: 

where the functions Cn{x,Q'^) depend weakly (i.e. logarithmically) on Q^. Various terms 
in this expansion are referred to as leading {n = 0) and higher {n > 1) twists. The 
"twist" number is defined in such a way that the leading one is equal to two and higher 
ones correspond to consecutive even integers. 



Thus the right hand side of the equation (|TT]) with approximate Bjorken scaling of 
quark distributions corresponds to the "leading twist" contribution to the F2. For mod- 
erately large values of (Q^ of the order of a few GeV^) contributions of the "higher 
twists" may become significant. Contrary to the common opinion the higher twists are 
only corrections to the leading (approximately scaling) term (0) in the large region. 
Thus they cannot correctly describe the low (i.e. nonperturbative) region since the 
expansion ([T2|) gives a divergent series there. In particular it should be noted that the 
individual terms in this expansion violate the constraint (^. In order to correctly describe 
this region the (formal) expansion has to be summed beforehand, at large Q^, and then 
continued to the region of ~ 0. This is automatically provided by certain models like 
the Vector Meson Dominance (VMD) model. To be precise the VMD model together 



4 



with its generalisation which gives (approximate) scahng at large can be represented 
in a form ([l^) for sufficiently large Q^. 

The VMD model is a quantitative realisation of the experimental fact that the photon 
interactions are often similar to these of a hadron [jlO|, |Tl|. The structure function F2 is 
in this model represented by: 



X 



QVis.Q^-M%Q^] ^ (13) 



The quantities (t^{s) are the vector meson-nucleon total cross sections, is the mass of 
the vector meson v and 7^ can be related in the standard way to the leptonic width of 



the V 10 



If only the finite number of vector mesons is included in the sum ( |T3D then the F2 vanishes 
as at large Q^. Therefore it does not contain the "leading twist" term. The scaling 
can be introduced by including the infinite number of vector mesons in the sum. This 
version of the VMD is called the Generalised Vector Meson Dominance (GVMD) model 



10, 11]. The heavy mesons contribution is directly related to the structure function in 



the scaling region. 

In practical applications to the analysis of experimental data which extend to the 
moderate values of one often includes the higher twists corrections in the following 
simplified way: 



1 + 



(15) 



where the F^'^ is the leading twist contribution to F2 and H{x) is determined from fit 
to the data. This simple minded expression may not be justified theoretically since in 
principle the higher twist terms, i.e. functions C„(x,(5^) for > 1 in eq.(12) evolve 
differently with Q"^ than the leading twist term. 



At high energies the Regge theory is often used to parametrise the cross sections 



as the functions of energy W . The high energy behaviour of the total cross sections is 
then given by the following expression: 

<^i{W) = Y.(3.{W^r-' (16) 

i 

where are the intercepts of the Regge poles and (3i denote their couplings. The intercepts 
Qfj are the universal quantities, i.e. they are independent of the external particles or 
currents and depend only on the quantum numbers of the Regge poles which are exchanged 
in the crossed channel. The Regge pole exchange is to a large extent a generalisation of 
the particle exchange and formally describes a pole of the partial wave amplitude in the 
crossed channel in the complex angular momentum plane. The Regge pole corresponding 
to the vacuum quantum numbers is called pomeron. It is expected to have the highest 
intercept, close to unity. It is a phenomenologically established fact that the energy 
dependence of the total hadronic and photoproduction cross sections can be described 
by two contributions: the (effective) pomeron with intercept Q;p=1.08 and the Reggeon 
with the intercept close to 0.5 [1^. Since bulk of the cross sections comes from the 
"soft" processes one usually refers to the phenomenological pomeron having its intercept 
slightly above unity as the "soft" pomeron. On the other hand in the leading logarithmic 
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approximation of perturbative QCD one finds the pomeron with intercept significantly 



above unity [|14|, |1^. Since perturbative QCD is only applicable for the description of the 
"hard" processes one refers to this QCD pomeron as the "hard" pomeron. Its relation to 
and interplay with the phenomenologically determined "soft" pomeron is still not fully 
understood. It should also be emphasised that the pomeron with intercept above unity 
leads to violation of unitarity at asymptotic energies. The unitarity is restored by multiple 
scattering absorptive (or shadowing) terms. 

The Regge parametrisation of the total cross sections implies the following parametri- 
sation of the electroproduction structure function F2(x,Q^): 

F,{x,Q^) = Y.m'){w'r-' (17) 

i 

It is expected to be valid in the high energy limit and for W'^ ^ Q^. In this limit 
X ~ Q'^/W'^ and x <^ 1. Therefore the Regge parametrisation (|T^ implies the following 
small X behaviour of F2{x, Q"^): 

F2(x,Q2) = EA(QV-"' (18) 

i 

where 

At low Q"^ the functions l3i{Q'^) should obey the constraint (||), i.e. l3i{Q'^) = Oi^Q"^) for 



4 High energy photoproduction 



The low region should join smoothly the photoproduction limit, = 0. In this limit 
the following relation between the total photoproduction cross section cr^p{E^) and the 
structure function F2 holds: 




This limit should be taken at fixed v = where is taken as photon energy in the 
laboratory frame. Several structure function parametrisations use the photoproduction 
total cross section as an additional constraint. 

When considering photoproduction it is conventional to decompose the total cross 
section ayp into two parts: 

0"7P = O'VMD + CTpart (21) 

and then 



^part ^direct ~l~ ^anomalous (22) 

In this equation avMD denotes the cross section corresponding to the VMD photon-hadron 
interaction mechanism while apart to the partonic mechanism respectively. The partonic 
component of the cross section is next decomposed into two terms: the "direct" term 
which reflects the photon interactions with the partonic constituents of the hadron and 
the "anomalous" term which corresponds to the interactions of the partons - constituents 
of the photon with partonic constituents of the hadron. In the latter case the photon 
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coupling to its constituents is point-like. The main feature of the partonic mechanism is 
that it corresponds to the (semi) hard interactions which can be described by perturbative 
QCD. The VMD part on the other hand contains both hard as well as the soft components. 
The former comes from the hard interaction of the partons which are the consituents 
of the vector meson. Unlike the "anomalous" part this term cannot be described by 
perturbative QCD alone since its description requires knowledge of the (nonperturbative) 
parton distributions in the vector meson. The anomalous term together with the hard part 
of the VMD contribution represent the point-like interaction of the partonic constituents 
of the photon. In the literature the corresponding events are therefore called the "resolved" 
photon events |16|, |1^, |18|, p^ ]. 



It follows from (pOD that the total photoproduction cross section can be obtained from 
the extrapolation of F2/Q^ to = 0. The partonic component <Jpart (i-e. the sum of 
the "direct" and the "anomalous" terms) can be identified with the difference between 
the result of the extrapolation and the VMD part |T^. The most important distinction 
between various mechanisms can be done through the analysis of the final states which 
will not be discussed here (see 

The photoproduction total cross section, like the hadron— hadron total cross sections 



exhibits an increase with the increasing photon energy, ^ see fig.3. This increase 
can be well described by the "soft" pomeron contribution with its intercept equal to 1.08 
1l3| (see also sec. 3). It is also well described by the extrapolation of the parametrisation 



by Abramowicz et al., |23]. Both predictions are shown in fig. 3. Another possible 



description of the total cross section increase is provided by the (mini)jet production, i.e. 
production of jets with relatively low pt (beginning from 1—2 GeV or so). The minijet 
production is described by the hard scattering of partons coming from the photon and 
from the proton respectively. The energy dependence of the cross section reflects to a 
large extent the small x behaviour of the parton (i.e. mostly gluon) distributions in a 
proton and in a photon. This follows from the kinematics of the hard parton-parton 
scattering: 

XiX2W^p > (23) 

where xi^2 denote the momentum fractions carried by partons and W^p is the total CM 
energy. For increasing W^p (and for fixed p^) , Xi and X2 may assume smaller values. 
The magnitude of the cross section is sensitive to the magnitude of the minimal value of 
Pt- It can be shown that the multiple scattering (or absorptive) corrections are extremely 
important here and they slow down significantly the increase of the total cross section with 



energy fl^, 0, Possible prediction for the total cross section energy dependence 
which follows from the minijet production picture are shown in fig. 3 (see |2^ for the 
details). 

The global quantity like the total cross section is not capable to discriminate between 
different models and the detailed structure of final states may be crucial here. The relevant 
Monte Carlo event generator for the minijet model with multiple scattering is discussed 
in [EH . 
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5 Phenomenological parametrisations of structure func- 
tions 



There exist several phenomenological parametrisations of the structure functions which 
incorporate the — constraints (cf. sec. 2) as well as the Bjorken scaling behaviour 
at large |]2^, |30|, |23|, |32|. Certain parametrisations ^ also contain 

the (QCD motivated) scaling violations. The parametric form of the corresponding 
dependence, however, is not constrained at large by the Altarelli-Parisi evolution 
equations, i.e. those parametrisations are not linked with the conventional QCD evolution. 
Nor is the low behaviour related to the explicit vector meson dominance, known to 
dominate at low Q^, [llO| . 



A parametrisation used by CHIO Collaboration |27| to fit their data is a com- 
bination of a simple parton model (valence quarks) and the GVMD spirit approach (sea 
quarks) : 

F2(x, Q2) = P3(2 + g,)x{\ - x)i+^^ + ^(1 - ^f^''7^^2 (24) 

9 5 + + mg 

where = go3 + e, g^ = go5 + e, e = Kln[{Q'^ + mg)/mg] and P3, P5, 5(03, fi-os, and mo 
are free parameters. This parametrisation is obsolete and has been mentioned here for 
historical reasons only. 

A parametrisation by Brasse et al. concerns the resonance region 0.1 < < 



6 GeV^, 1.11 < W <1.99 GeV, cf. fig.l. The virtual Compton scattering cross section is 
assumed as: 

a(g^ W, e) = ariQ^ W) + ea^Q', W) (25) 

where the parameter e in the formula (^) denotes the degree of polarisation of the virtual 
photon and the cross sections ariQ"^, W) and (Ji(Q^, W) are related in a standard way to 
the structure functions F2 and Fj^: 



ctt + ctl) (26) 
ai (27) 



The cross section cr(Q^, W, e) is parametrised in the following way: 

ln{a/G^) = a{W) + b{W)ln}-^ + c{W) \ In^^ l'^(^) (28) 

where q is the three momentum transfer to the hadronic system , i.e. | |= \/Q^ + i^'^, 
I Qo I is the value of | q | for Q'^=0 , i.e. | Qo \= (W^ - M^)/2M and G^{Q^) is the dipole 
form factor of the nucleon, i.e.: 

The parameters of a{W), b{W), c{W) and d(W) were obtained from the fit to the data in 
different bins of W. Their values are tabulated in ref . ||28[] . 

A complete parametrisation of P^ including the resonance region was ob- 
tained by the NMC [^]. In the resonance region F2 was fitted to the data from SLAC 
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P^ l taking only the A(1232) resonance into account. Outside the resonance region a QCD 
based parametrisation was used to describe the data of OHIO SLAG @], BCDMS 
p5| and EMC NA28 For this purpose, the structure function was parametrised as: 



(30) 



where F'^*'^ and F^*^* are the contributions from the deep inelastic and resonance regions 
respectively and describes the background under the resonance. The nucleon electro- 
magnetic form factor is given by eq. (^9]) ; the factor 1— in eq. (|30| ) suppresses F2 at 
low values of where elastic scattering on the nucleon dominates. 



The contribution from the deep inelastic region was parametrised as 



:^d^s , 



18 B{r],,r]2 + 1 



<(1 



VI4 



Six,Q') (31) 



where x^ = {Q^ + ml)/{2mu + ml) with m^=0.351 GeV^ and ml =1.512 GeV^. The 
quantity B is the Euler's beta function and 771, 774 are linear functions of the variable s, 



Oii + A-s, 



(32) 



where 



In 



ln[{Q^ + ml)/A^ 



(33) 

04 and /?4 were free 



'ln[{Ql + ml)/K^] 

with Qg=2.0 GeV^ and A=0.2 GeV. The constants ai, 
parameters in the fit. 

The factor ^(x, Q^) in eq.(|3TD suppresses F'^^^ in the resonance region close to the 
single pion production threshold: 



^(x,Q2) = l-e-''(^-^'-), 

with Wthr =1.03 GeV and a=4.177 GeV'^ 

The form adopted for the contribution from the resonance region was 



(34) 



(35) 



with mA=l-232 GeV, r=0.0728 GeV and 05, a free parameter in the fit, was equal to 
0.89456. This parametrisation takes into account only the A(1232) contribution; higher 
mass resonances are neglected. 



The background under the resonance region was parametrised as 



F'^{x,Q') = alG'"ie 



■b{W-Wthr? 



(36) 



where 
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\ 



[{W + cY + M"^ - mlf 



M2 



(37) 



with b=0.5 GeV ^ and c=0.05 GeV. The parameter ag, left free in the fit, was equal to 
0.16452. 
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F2{x,Q^) = A{x) 



2 1 \2\lB{x) 



ln(gVA 



1 + 



C{x) 



Qg = 20GeV^ A = 250MeV; 
A{x) = (1 - x)"^ [as + 04(1 -x) + 05(1 - xy + 06(1 - xf + aril - x^] ; 

B{x) = bi + b2X + bs/ix + 64); 
C{x) = CiX + C2X^ + c^x^ + c^x'^. 



Parameter 


proton 


deuteron 


ai 


-0.1011 


-0.0996 


0-2 


2.562 


2.489 


03 


0.4121 


0.4684 


04 


-0.518 


-1.924 


as 


5.967 


8.159 


oe 


-10.197 


-10.893 


aj 


4.685 


4.535 


bi 


0.364 


0.252 


b2 


-2.764 


-2.713 


bs 


0.0150 


0.0254 


&4 


0.0186 


0.0299 


Cl 


-1.179 


-1.221 


C2 


8.24 


7.50 


C3 


-36.36 


-30.49 


C4 


47.76 


40.23 



Table 1: The parametrisation of F2 and j^. This function is strictly valid only in the 
kinematic range of the NMC, SLAC and BCDMS data. 



The parametrization (^) of the function F2 is vahd from Q ~ up to about 200 



GeV^ and from x = 0.003 up to 0.7. However, the results of the fit given in 



are 



obsolete since the new F2 measurements by the NMC [Q which were not included in the 
fit differ by up to 30% from the former world data. When used later in the radiative 
corrections procedure the deep inelastic part of F2 was refitted while the parametrisation 
of the resonances and of the background was kept fixed. Two types of the former were 
used: 8- and 15-parameter functions. The 15-parameter function and the fitted values 
of its coefficients are given in Table 1 [0]. 

A low parametrisation of was obtained by Donnachie and LandshofF 
PD[ . In this parametrisation the parton model is extrapolated down to the region of 



low Q"^ (including the photoproduction) respecting the constraints discussed in sec. 2. 
The parametrisations of the sea- and valence quark distributions at small values of x 
are based on Regge theory (cf. sec. 3) and correspond to the pomeron intercept equal 
approximately 1.08 and reggeon intercept about 0.5. The structure function is given 
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by the conventional parton model formula: 



= x[^{u + u) + ^{d + d) + ^X{s + s)] (38) 

where u, u,...aTe the parton densities ans A measures the reduced strength of the pomeron's 
couphng to strange quarks (A ~0.6). The individual parton densities are related in the 
following way to the sea- and valence quark distributions, S{x) and V{x) respectively. 

xu{x) = S{x) + 2V{x), xd{x) = S{x) + V{x), xs{x) = S{x) 

xu{x) = xd{x) = xs{x) = S{x) (39) 

The densities S{x) and V{x) are parametrised as follows: 

S{x) = 0.17a;-°■°^(l-a;)^ 

V{x) = 1.33x°-^^(l -x)^ (40) 

The functions S{x) and V{x) are next multiplied by the factors 0s(Q^) and (p^iQ"^), 
respectively, where 



+ 0.36 

' Q 



2 \ 0-44 



g2 + 0.85 



(41) 

The powers are chosen so as to make vanish like as 0. This parametrisation 

has been recently extended and modified ||31| to include in more detail the heavy quark 
contribution and the region of large x. Instead of the scaling variable x one uses: 

f.=.(l + |) (42) 

where the index 'i' denotes charm, strange and light quarks. Strange and charm distri- 
butions are as follows: 

xs{x,Q^) = c.-^er(i-es)' 

xc{x,Q') = C,-^er(l-ec)' (43) 

where ^ 0.22, =0.032, a, =1 GeV^, = 6.25 GeV^ and e = 0.0808. Parameters 
^5^^ defining the variables C,s,c (cf. eq. (|4^) ) are: fi^ =1.7 GeV^ and fi^ =16 GeV^. 

Parametrisation of the light quark distributions is more involved, i.e. different func- 
tional forms obtained from a fit to the data are used for C, smaller and greater than 
a certain which is a free parameter in the fit. The valence quark distributions are 
parametrised as below: 

xu,ix,Q'') = f/(OV'(Q') 

UiO = B^e'''' C<Co (44) 
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and 

^(0 = Pde^a-o' e>eo (45) 

The parameter /i^ defining the variable ^ for hght quarks is /i^ =0.28 GeV^. Parameters 
Bu,d, Pu,d and Xu,d are fixed by requiring that the two functional forms defining U{^) 
and D{C,) join smoothly at and by imposing the number sum rules. The function ip is 
i^iQ"^) = Q'^/iQ'^ + b)- The light sea quark distributions are: 



qsix,Q') = ^e-il-O'm') e>eo (46) 

where $(Q^) = Q"^/ {Q"^ + a). The parameters C, a and b are fixed by fitting to the data 
on photo- and electroproduction. The model contains also a Regge-like term in the sea 
quark distributions which behaves like ^"^-^^^s small x. It also contains an extra "higher 
twist" term in the F2 which is parametrised as below: 



The values of other parameters defining the model are given in The parametrisation 
is vahd for 0< <10 GeV^ and for x >0.008. 

The function -^^'(x, Q^) resulting from the parametrisation of Donnachie and Landshoff 
is shown in fig. 4 together with the data. 



A parametrisation proposed by Abramowicz et al., 23 1, is a result of a fit to 



the following structure function data sets: SLAG H, BGDMS H and EMC NA28 H 
and the then available photoproduction cross section. The fitted function is based on the 
parton model with the QGD motivated scaling violation appropriately extrapolated to 
the region of low (including the photoproduction). The Regge ideas are used in the 
parametrisation . 

The structure function F2 is decomposed into two terms, F2 and F^, corresponding 
to pomeron and reggeon exchange, 

F2 = Fi + Fi- (48) 
with each of them expressed as follows: 



P2 = „ , a(t)x^-(*)(l - x)'^(*) (49) 



where 

1 2Mi/ + m2 , ^ 

Xr + 

and r = P,R stands either for pomeron or for reggeon. The argument of the functions 
Crit),ar{t) and br(t) is defined as: 
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The functions Cr, bji, and hp which increase with Q"^ were assumed to be of the form: 



pit) =Vi + {Vi-V2y' (52) 
while the functions Cp and ap were assumed in the form: 

pit) =pi + (pi - P2) - 1) (53) 

The parameters mQ,m'^,Ql as well as the parameters ai,bi and Ci (separately for the 
pomeron and the reggeon terms, 2=1,2,3) were obtained from the fit to the data. 

The fit of Abramowicz et al. is valid from =0 up to the highest values obtained 
in the fixed target experiments, excluding the resonance region, W <1.75 GeV. The fit 



given in |^ should, however, be updated since the recent NMC data on F2 [|1| were 
not included. The structure function (x, Q^) which follows from this representation is 
shown in fig. 5. 



A parametrisation by Capella et al. is based on a fit to the fixed target data as 



well as to the most recent large data from HERA, P7| , 3H| which show the significant 
increase of F2 with decreasing x. It also uses the data on the photoproduction cross 
section including the highest energy results from HERA [^, |2^. This parametrisation 



assumes a Regge form of F2 at small x and counting rules at large x. The leading small 
X behaviour corresponding to the pomeron is parametrised by the dependent effective 
pomeron intercept 1+A((5^) which interpolates between the "soft" pomeron at Q^=0 and 
the "hard" one at large Q^. The unusual dependence of the effective pomeron intercept 
is meant to describe in a phenomenological way the absorptive corrections to F2 which 
are expected to be strong at low values. The small x behaviour corresponding to the 
Reggeon contribution is parametrised in a conventional way, cf.eqs (p!8|JT9| ). To be precise 
the structure function F2{x,Q'^) is parametrised in the following form: 

F2{x,Q') = Ax-^(«^)(l-x)-(«^)+M-^ 

+ a) 

+ Sa;i-"«(l-x)"(«'M (54) 

The first term in (^4]) corresponds to the pomeron contribution where the effective pomeron 
intercept is parametrised as: 

with Ao ~ 0.08 that corresponds to the intercept of the "soft" pomeron. The second term 
in (|5^) corresponds to the reggeon contribution and at large can be identified with the 
valence quarks. The pomeron part corresponds to the sea quark contribution at large . 
The function n{Q'^) is parametrised as follows: 

Parameters A,B,a, b, c, d are obtained from the fit to the data. The parametrisation 
(|5^ is used for moderately small and small (including photoproduction, i.e. Q^=0). 
The region of large is described by the combination of this parametrisation with the 
Altarelli-Parisi evolution equations. 
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6 Dynamical models of the low behaviour of 



In the previous section we have presented various parametrisations of the structure func- 
tions. Those parametrisations although being motivated dynamically were not linked 
with conventional QCD evolution at large (except of that by Capella et al., |32]). 
Nor was the low dynamics explicitely taken into account. In this section we present 
parametrisations which will contain the QCD evolution. One of them will also include 
the VMD dynamics, dominating at low . 

Dynamically calculated structure functions by Gliick, Reya and Vogt ||39|, 
extend the QCD evolution equations down to the very low region (Q^ < 1 
GeV^). The parton distributions are calculated by evolving from the input valence-like 
distributions provided at the very low scale /x^ =0.3 GeV^ and both the leading and 
next-to-leading order approximations are used. The valence-like shape at the scale ^ is 
assumed not only for the valence quarks but also for gluon and sea quark distributions: 



xgKx,[i-) = Ax''{\-xf , xq{x,ii') = Nx"^ {l-xf (57) 

The parametrisation of the valence quarks is provided at the large scale =10 GeV^ 
and its form at = is obtained from the backward QCD evolution. Parameters 
describing the parton distributions come from a fit to the data at large and are also 
constrained by the sum rules. As a result of the QCD evolution the gluon and sea quark 
distributions (multiplied by x) become immediately singular at small x for > fi"^. In 
the leading twist approximation the absolute value of the slope in x of these distributions 
and so the absolute value of the slope in x of the F2 grows with increasing 

Q2. Naturally 

the model which is based on perturbative QCD cannot be extended into the region of very 
low < where A is the QCD scale parameter since the (perturbatively calculated) 
QCD coupling as{Q'^) 00 for A^. Moreover the measurable quantities like the F2 

may contain non-negligible higher twist contribution in the region of moderately small 
Q^, cf. eq. ([T^). Presence of higher twists in an observable should not affect the QCD 
evolution of the leading twist parton distribution since partonic distributions of different 
twists evolve separately. Surprisingly the QCD evolution turns out to be stable and gives 
positive definite parton distributions down to the very small scale = (2A)^. Also it 
turns out that in the next-to-leading order approximation of perturbative QCD, the F2 is 
more stable than the distributions of partons |41|. Although the higher twist contribution 
can in principle be important in the low region, the leading twist F2 calculated within 
the model has recently been succesfuUy confronted |^0[ with the low x and moderate 
(Q^ equal few GeV^) NMC data. The high Q^, low x results of HERA are resonably 
described too. 

The explicit forms of the parton distributions' parametrisations are as follows: 

xv{x,Q^) = Nx''{l + A^/^+ Bx){l-x)^ (58) 

for valence quarks 



A'x"'(l 



xw{x, 



x"(A + Bx + Cx^) ( In^ + s^exp (-E + J E'sH'n}- 



for gluon and light sea quarks and 

is - sS"' 



{l-xf (59) 



xw'{x, Q 



1 + ^Bx)(\-xY exp \-E + d E'sf^ln 



X 



(60) 
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for heavy sea quarks, where s in the formulae (BO-pO) is: 



In these equations A, B, C, D, E, E', a and b are low order polynomials of s. These poly- 
nomials as well as values of remaining parameters can be found in . 



Electroproduction structure function F2 in the low Q , low x region by 



Badelek and Kwieciiiski p2| , Contributions from both the parton model with 
QCD corrections suitably extended to the low region and from the low mass vector 
mesons were taken into account. The former contributions results from the large 
structure function analysis [0, ^ which includes the recent F2 measurements by the 
NMC 0. 

The starting point is the Generalised Vector Meson Dominance (GVMD) representa- 
tion of the structure function F2{x,Q'^) [Q: 

^ Ft\x,Q') + Fi'\x,Q') (62) 
The function $((5^, s) is expressed as follows: 



$(Q^s) = --/m/ - :^F,^'{x',Q") (63) 



1 ^ f~Q' dQ'' 
—Im / — 
TT J Q'2 

The asymptotic structure function F^^{x,Q'^) is assumed to be given. It may be ob- 
tained from the QCD structure function analysis in the large region. By construction, 
F2{x,Q^) — i> F^^{x,Q'^) for large Q^. The second term in (p^) can be looked upon as 
the extrapolation of the (QCD improved) parton model for arbitrary Q^. The first term 
corresponds to the low mass vector meson dominance part since the sum extends over the 
low mass vector mesons. Contribution of vector mesons heavier then Qq is included in 
the integral in (|^). Choosing the parameter Ql > {M^)max where {Mv)max is the mass 
of the heaviest vector meson included in the sum one explicitly avoids double counting 
when adding two separate contributions to the structure function. Note that Qq should 
be smaller than the mass of the lightest vector meson not included in the sum. The 
representation (|6^ is written for fixed s and is expected to be valid at s Q^, i.e. at 
low X but for arbitrary Q^. 

In the representation ( |62D for the partonic part F2'\x,Q'^) was simplified as 
follows: 

F^'\x, Q') = jQ^^/2^'i^^ Q" + ^0) (64) 

where 

Q' + Ql _ Q' + Ql 

Simplified parametrisation ( |64D connecting F2^\x,Q'^) to F.^^ by an appropriate change 
of the arguments of the latter posseses all the main properties of the second term in 
(H). First of all it is evident that F^^\x,Q'^) F^^{x,Q^) for large Q"^. Moreover the 
parametrisation of F^"^ defined by (|6^ ) preserves the analytic properties of the second 
term in the eq.(B^). 
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It should be stressed that apart from the parameter Qq which is constrained by physical 
requirements described above the representation (^) does not contain any other free 
parameters except of course those which are implicitly present in the parametrisation of 
parton distributions defining ^2'^'^. 

The proton structure function calculated from the representation (|62|) with ^2''^ [x, Q^) 
given by the eq.(64) is shown in fig. 6. In an "approach to scaling" (i.e. in the de- 
pendence of F2 at less than 1 GeV^ or so) visible in fig. 6 the change of curvature 



comes from the factor Q'^ in eq. (|6^) and the magnitude of this change, particularity at 
^ 1 GeV^, is controlled by the vector meson contribution, a non-trivial test of this 
mechanism. Although the vector meson contribution Fj^"* dominates in the very low 
region (Q^ < 1 GeV^) the partonic component Fg^"* still gives a significant contribution 
(i.e. at least 20% - 30%) there. At = 10 GeV^ the structure function F2 calculated 
from the model differs from F^^{x, Q^) by less than 3% . 

The X dependence at fixed reflects both the energy dependence of and the 
a:-dependence of F^^. Expectations coming from Regge theory are incorporated in the 
parametrisations of the a^^s), and of the input parton distributions at the reference 
scale = 4 GeV^, 1^^. F2{x,Q'^) increases with decreasing x and that reflects the 
increase of the total cross sections cr^,(s) with increasing s as well as the increase of 
F2'\x,Q'^) with decreasing x. The increase of F2'\x,Q'^) with decreasing x is weaker at 
small than the similar increase in the large region. The s dependence of ay{s) is 
also relatively weak in the relevant region of s, As the result the x dependence of 
F2{x, Q"^) is relatively weak for low Q^. 

The model is applicable in the small x region {x <0.1) and for arbitrary value of Q^, 
however the way of parametrising the total vector meson-nucleon cross sections imposes 
additional constraint: u >10 GeV. 



A parametrisation by Schuler and Sjostrand extends those by Abramowicz 
et al. and by Donnachie and Landshoff including explicitly the damping of structure 
functions at small x which is implied by screening effects in QCD JT3| iB| and taking 



into account the logarithmic scaling violations at large Q^. Extrapolations towards the 
region of small are done applying the same damping factors as in refs. |^ |3^, 
In the region of low x the Regge parametrisation is used. In practice the (x, /i^) plane 
(yU^ = Q^) is divided into four regions as shown in fig. 7 and different parametrisation of 
parton distributions is assumed separately in each region. For large /i^ (/x^ > /Zg) the 
regions of large and small x (regions I and III in fig. 7) are divided by the boundary curve 
/i|(x) 

/i|(x) = 2 + 0.053^exp ( 3.56 J In— 1 . (66) 




In the region I of large fi^ and large x the structure functions F2 is given by the QCD 
improved parton model with the scale dependent quark distributions fv,s{x, /i^) where the 
subscripts v and s denote the valence and sea distributions. The distributions fv,s{.x,jj?) 
satisfy the Altarelli-Parisi evolution equations. In the regions II-IV the structure function 
F2 is again expressed in terms of the quark distributions /„,s(x, /i^) with their /i^ depen- 
dence being different in different regions separately. In the region II of small /i^ and large 
X [x > xq) the structure functions are damped by the /i^ dependent factors and valence 
distributions are made harder, i.e. 

/ ,,2 ,,2 I ™2 \ {l-r,){^-x)/{l'Xo) 
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fsix, /Xq) 



(6^ 



In the region III of large /i^(/io < /i^ < fi%{x)) and small x, the sea and valence quark 
distributions are parametrised as below: 



xfs{x, /i^) = ( — ) Xofs{xo, l4) + N2Xfs{x, 

^Xq , 



(69) 
(70) 



where 



(71) 
(72) 



ln{fil/fil) 
N2 = l-Ni 

Finally in the region IV of small x and small yU^ the valence and sea quark distributions are 
parametrised in the Regge form with the appropriate damping factors which guarantee 
vanishing of the structure function F2 in the limit /i^ = 



Xfy{X,IJ,'^) = — 



X \i 

Xo. 



(73) 



X fs {X-i /i 



X 



iVi (^^±^\ xofsixo, fil) + N2X^'N' 
\ f^o J 



where now 



A^i = 1 - 



No 



2 2 



Xo 



(74) 



(75) 
(76) 



The parameters e and 77 correspond to the (effective) intercepts of the Pomeron and 
the Reggeon respectively (e = 0.56, rj = 0.45). The values of the remaining parameters 
appearing in the formulae ( |67D - (|75|) are as follows: /ig = 5 GeV^, Xq = 0.0069, mfj = 0.92 
GeV^,ml = 0.38 GeV^,N^ = 0.121, A^^ = 2iV^, iVJ = A^^ = = = 0.044 and 
= N? = 0.5N^. The structure function F2{x,Q'^) is expressed in terms of the quark 
(and antiquark) effective distribution functions fi{x, Q"^) using the parton model formula: 



F2(x,Q2)=x^e^/,(x,Q2). 



(77) 



Comparison of different structure function parametrisations is shown in fig. 8. 



7 Nuclear shadowing 

Nuclear effects in the deep inelastic lepton scattering may be experimentally inferred from 
inspecting the ratio F2/F2. In this method one tacitly assumes that the nuclear effects 
in F2 can be neglected or that = F2 where F2 describes a free nucleon. In particular 
the nuclear shadowing corresponds to the F2 / F2 ratio being smaller than unity at small 
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X. In this sense shadowing can be regarded as a part of the so called EMC effect as 
illustrated in fig.9. 

The nuclear shadowing is a firmly established experimental fact, cf.sec.8, which has 
been observed both in the low (including the photoproduction, i.e. = 0) and in 
the large region. The nuclear structure function is then: 

F,^ = F^ - 8F^ (78) 

where 5F2 denotes the shadowing term {5F2 > 0). 

In the region of low (and for photoproduction) the natural and presumably the 
dominant mechanism of nuclear shadowing is the multiple scattering of vector mesons in 
the nucleus, fig. 10. The vector mesons couple to virtual (or real) photons. This model 
gives the following contribution of the shadowing to the nuclear structure function: 

where the cross section 6a^ is that part of the vector meson - nucleus total cross section 
(normalised to a nucleon) which corresponds to multiple scattering i.e.: 

= - S^t- (80) 

Here is the vector meson - nucleon total cross section and the remaining quantities 
in the formula (I79|) are defined in sec. 3. The cross section which corresponds to the 



multiple scattering can be obtained from the Glauber theory [^, 0. The negative sign 
in eq.(pUD refiects the fact that the vector meson-nucleon scattering amplitude is assumed 
to be imaginary. Thus the shadowing in the inelastic lepton-nucleus scattering, in the 
VMD refiects the absorptive character of the elementary vector meson-nucleon scattering 
amplitude. It follows from eq.(^) that the shadowing term which corresponds to the 
multiple rescattering of (finite number of) vector mesons vanishes for large Q^. 

At large one expects the parton model to be applicable. It leads to the Bjorken scal- 
ing mildly violated by the perturbative QCD corrections. The parton model is described 
by the "hand-bag" diagram of fig. 11. It is this "hand-bag" structure and the point-like 
coupling of the photons to partons (i.e. to quarks and antiquarks) which guarantees the 
Bjorken scaling (modulo the perturbative QCD corrections) independently of the struc- 
ture of the lower part of the diagram. The nuclear shadowing in the large region may 
come from the multiple interaction contributions to the lower part of the hand-bag dia- 
gram as shown in fig. 11. Different models of shadowing correspond to different structure 
details of the diagrams of fig. 11 (see ^9]). 

Shadowing is expected to be a low x phenomenon. This can be understood within the 
simple space-time picture of the interaction of the virtual photon with atomic nucleus. 



pD| |. In the infinite momentum frame i.e. in a frame where the momentum p = Pa/A is 
very large the wee partons (sea quarks and gluons) occupy longitudinal distances 

Azp ^ — (81) 
xp 

The momentum pA is the momentum of the nucleus with A nucleons. The nucleus in this 
frame occupies the Lorentz contracted distance: 

M 

Aza ~ 2Ra— (82) 
p 
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where Ra is the nuclear radius and M the nucleon mass. The (Lorentz contracted) average 
distance between the nucleons in this frame is: 

M 

Az = r— (83) 
P 

where r denotes the average distance between nucleons in the laboratory frame. One can 
distinguish three regions in x, i.e.: 

(84) 

that corresponds to the partonic size being smaller than the average distance between 
nucleons within nuclei: 

Azp < Az. (85) 

In this region the shadowing is expected to be negligible since partons can be regarded as 
belonging to individual nucleons. 

(^^) 2MRa ^ ^ 

(86) 

that corresponds to the longitudinal size of the partons being larger than avarge distance 
between nucleons in a nucleus yet smaller than the (contracted) longitudinal size of the 
nucleus: 

Aza > Azp > Az. (87) 
In this region the shadowing gradually sets on with the decreasing x. 

(^^^) ^ < 2MRa 

(88) 

that corresponds to 

Azp > Aza (89) 

where shadowing is expected to be maximal. In the regions (i) and (u) the partons can no 
longer be regarded as belonging to individual nucleons since their longitudinal distances 
are larger than the average distances between nucleons exceeding eventually the size of 
the nucleus. 

Let us also notice that when the deep inelastic scattering is considered in the laboratory 
frame then the three regions of x, i.e. [i) — [Hi) are defined through the mutual relations 
between the lifetime of the qq fluctuation of the virtual photon: 

r-. = (90) 

and the characteristic distances r and Ra- The nuclear shadowing indeed is a typical 
small X phenomenon since l/(Mr) 0.1. 

There exist several models of shadowing incorporating the VMD and/or the partonic 
mechanisms. For a complete review of these models and a comprehensive list of refer- 
ences see ref. [|9||. Since that article has been completed new papers concerning the 



shadowing in the deuteron have appeared several of them inspired by the 

new measurements of the F2/F2 ratio (cf. sec. 8). Understanding of shadowing effects in 
the deuteron has become relevant in view of the increased precision of the measurements 
of the F2 / F2 and F^ / F2 ratios. In particular the shadowing affects extraction of the 
neutron structure function from the data. This in turn affects (decreases) the magnitude 
of the experimentally evaluated Gottfried sum ([Rl, |51[]). 
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8 Experimental data 



The experimental problems connected to measuring and analysing the low x, low data 
are discussed in ||^. In that reference we have also listed the experiments which provided 
the then available experimental results. Here we limit ourselves only to updating this 
information. 

In the fixed target experiments the low x region is correlated with the low values, 



see e.g. fig. 12, |5^. The lowest values of x were reached by the NMC at CERN and E665 
Collaboration at Fermilab through applying special experimental techniques permitting 
measurements of muon scattering angles as low as 1 mrad, [^6|, |5^. These "small 

X triggers" and special off-line selection methods were also effective against a background 
of muons scattered elastically from target atomic electrons; corresponding peak occurs 
at X =0.000545. Systematic errors in both experiments, in particular these on the ratio 
of structure functions for different nuclei, F2/F2, were greatly reduced as a result of 
irradiating several target materials at a time and/or of a frequent exchange of targets in 
the beam. 

During the last three years there appeared new measurements of the -Ff Q^); -^2^(^; Q"^) 
and F^{x,Q^)/F^{x,Q^) [0, |9|, ^ |3 by the NMC and by the E665 [0, ||]. The 
NMC also performed the QCD analysis of their F2 data 0. Both collaborations have 
presented new and precise results on x, A and dependence of nuclear shadowing 
|6|, H, The NMC also determined differences R^" - , and R"^ - Rp, || 

at low X. 

Extraction of F2{x, Q^) from the data needs information on R{x, Q"^). In particular the 
ratio of inelastic cross sections on different nuclei is equal to the corresponding structure 
functions ratio, provided R{x, Q"^) is the same for these nuclei. Results of the NMC 
analysis of i?*-^" — i?*^ and R'^ — R^ show that neither of these quantities exhibit a significant 
dependence on x and that they are both compatible with zero, fig. 13. In their analyses 
the NMC and E665 assumed R independent of the target atomic mass A. 

The NMC results for the deuteron structure function, as a function of for 
different bins of x, are shown in fig. 14. Measurements cover ranges 0.006 < x < 0.6 and 
0.5 < < 55 GeV^ and were the first precise measurements at such low values of x. The 
data had great impact on the parton distribution analysis (see e.g.|^^) and join well to the 



results of HERA, A clear scaling violation pattern with slopes d lnF2/d InQ^ positive 
at low X and negative at higher x and an "approach to scaling" (i.e. dependence of 
F2 at less than few GeV^) are visible in fig. 14. In this figure comparison of the NMC 



[]1|, SLAC and BCDMS measurements is also shown. The agreement between 
all three data sets is very good. At the same time the low x results of the EMC NA2 
p7 | experiment have been disproved by the NMC measurements. The x dependence of 



the deuteron structure function at low is shown in fig. 15 for several values of Q^. 
The NMC measurements are compared with those of EMC NA28 and SLAC |^ . 
Characteristic is a weak x dependence of F2 at low x. 

The dependence of the structure functions and F2 measured by the NMC 
with good accuracy down to low values of x has been compared with the predictions of 
perturbative QCD [Q. The flavour singlet and non-singlet quark distributions as well as 
the gluon distribution have been parametrised at the reference scale equal to 7 GeV^. All 
the data with > 1 GeV^ were included in the fit. Besides the leading twist contribution 
the higher twist term was also included in an approximate way given by formula ( P^ ) where 
H{x) was determined from the SLAC and the BCDMS measurements averaged over 
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the proton and deuteron and suitably extrapolated to lower values of x. Results of the 
QCD fit to the proton structure function data are shown in fig. 16. Important here is 
the extension of the QCD analysis to the low x and low regions. The contribution of 
higher twists is still moderate at scales about 1 GeV^. 

The new (preliminary) measurements of the proton and deuteron structure functions 



for X >0.0001 have recently been presented by the E665 Collaboration, fig.17 |58|. The 
lowest values in their data reaches a very small value of 0.2 GeV^. 

Both NMC and E665 experiments have measured the deuteron to proton structure 
function ratio, F^/Ff , extending down to very low values of x. In case of the NMC the 
ratio has been measured directly, i.e. the measurement of the absolute structure function 
is used only for the radiative corrections calculations. The data are usually presented as 
the ratio jF^ where Fll is defined as 2F^ — Fj* . This quantity would give the structure 
function of the free nucleon in the absence of nuclear effects in the deuteron. NMC results 
for F2IF2 as the function of x are shown in fig. 18, E665 results are presented in 
fig. 19, [^. In the latter case instead of the ratio Fl'/Ff the ratio an jap is given. It 



is equal to the former due to A— independence of R. In both data sets the average 
varies from bin to bin reaching down to < >= 0.2 GeV^ at x =0.0008 for the NMC 
and < >= 0.004 GeV^ at x = 5x 10^^ for E665. The results of both experiments 
show that the F2/F2 stays always below unity down to the smallest measured values of 



X which at low x can be attributed to the nuclear shadowing in the deuteron [51, 52, 53 



New data appeared on the nuclear shadowing. Preliminary results of a high precision 
study of the A dependence of nuclear shadowing by the NMC performed in the range 
0.01 < X < 0.7 and 2 GeV^ < < QO GeV^ are shown in figures 20-22. The structure 
function ratios F2 /F2 for A = Be, Al, Fe and Sn together with earlier data of SLAG p9| , 
show a detailed pattern of the x dependence of shadowing, fig. 20. The NMC data cover the 
A range from A = 2to A = 118. In fig. 21 they are shown as a function of A for three bins 
of X. The functional dependence of F2 / F2 on A is approximately logarithmic (except, 
possibly, the light nuclei) and the structure function ratio has therefore been parametrised 
as F2 / F2 = a + blnA in each bin of x. The slopes, b, are displayed in fig. 22. The amount 
of shadowing increases strongly with the mass number A. Much lower values of x and 

are covered by the nuclear data obtained by the E665: x >0.0001 and >0.1 GeV^ 
thanks to a special trigger and off-line analyses, |6^, fig. 23. Shadowing seems to saturate 



at X about 0.004 as was observed earlier by the E665 [31] and indicated by the preliminary 
NMC data on the F^'/F^ and F^/F^ ratios measured down to x =0.0001 and (5^=0.03 
GeV^, [|55|. No clear dependence is visible in the E665 data in a wide interval of Q^, 
fig.24. 



9 Conclusions and outlook 

In this review paper we have summarised the present understanding of the electropro- 
duction structure functions in the region of low values of x and Q^, which have recently 
been measured in the fixed target, charged lepton inelastic scattering experiments. This 
has included a survey of theoretical constraints and expectations, clarification of certain 
definitions and concepts, collection of the existing parametrisations of structure functions 
in this region and presentation of the experimental data. We consider this paper as being 
an extension of the comprehensive review article on the small x physics [0| towards the 
more detailed treatment of the low problems. This means in particular that in the 
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summary of the experimental data we have presented only the results which appeared 
after the above mentioned article had been completed. 

Important property of the structure function F2{x,Q'^) which follows from the con- 
servation of the electromagnetic current is its linear vanishing as the function of for 
(for fixed u). This means that Bjorken scaling cannot be a valid concept at low 

and so the pure partonic description of inelastic lepton scattering has to break down 
for moderate and low values. At moderately large the higher twist contributions to 
F2 which vanish as negative powers of are frequently being included in the QCD data 
analysis. One also expects that at low the VMD mechanism should play an important 
role. 

The small x behaviour of the structure function F2{x, Q"^) is dominated by the pomeron 
exchange. Analysis of the structure function in the small x region for both low and mod- 
erate values of can clarify our understanding of the pomeron, i.e. possible interplay 
between the "soft" and "hard" pomerons, the role of shadowing (or absorptive) correc- 
tions, etc. At large the problem is linked with the QCD expectations concerning the 



deep inelastic scattering at small x, see ref.|70[ for a recent review. Besides the structure 



functions (or total cross sections) complementary information on the pomeron can also 
be obtained from the analysis of diffractive processes in the electro- and photoproduction. 
This concerns both the inclusive diffractive processes and the diffractive production of 
vector mesons; the present experimental situation is described in refs. ^ 

Descriptions of the low Q^, low x behaviour of F2 range from pure fits to experi- 
mental data to the dynamically motivated models. The existing parametrisations have 
been collected in a way which should make it easy to use in practical applications, e.g. 
in the radiative correction procedure. Wherever possible a dynamical content of each 
parametrisation was exposed. 

Since 1992 a wealth of measurements of F2 has been published. In the region of interest 
for this paper this included the NMC and E665 results extending down to very low x 
and values and displaying characteristic "approach to scaling" behaviour. Nuclear 
shadowing phenomenon was studied in great detail for targets ranging from A = 2 to 
A = 118 by the same two collaborations. Its x, and A dependence was precisely 
measured. 

The new possibilities have opened up with the advent of the HERA collider. Presently 
the data there are collected at larger than approximately 5 GeV^ and therefore have 
not been discussed in this review. These data show a very strong increase of F2 with 
decreasing x, p^, p8| . Let us remind that for fixed the increase of the structure function 



F2{x,Q'^) with decreasing x reflects the increase of the virtual Compton scattering total 
cross section with W"^. We would like to emphasise that this increase observed at large 
is much stronger than the increase of the total photoproduction cross section observed 
between the fixed target and the HERA energies. Dynamical understanding of this evident 
variation of the dependence with is certainly very interesting and important. The 
experimental data at low values of {Q^ ~ 1 GeV^) which would cover similar range of 
W"^ as those at high Q"^ would be extremely valuable for this purpose. It would therefore 
be very interesting to extend the HERA measurement down to the low values of as it 
has been planned recently |]73 . 
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Figure Captions 

1. Illustration of the continuity of physical processes: double differential cross section 
for the electron-proton inelastic scattering is sketched as a function of the energy 
transfer u for different values of the resolution Q^. Dashed and continuous lines 
correspond to constant values of x and W respectively. Definitions of kinematical 
variables are given in sec. 2. 

2. Kinematics of inelastic charged lepton-proton scattering in the one photon exchange 
approximation and its relation through the optical theorem to Compton scattering 
for the virtual photon. 
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Total photoproduction cross section as a function of the 7p centre of mass energy, 
W^p. Low energy data come from ref.[^, recent ZEUS Collaboration data - from 



ref. [^], previous data - from ref.[^. The curves are explained in the text. Figure 



comes from ref. 22 



4. Parametrisation of F2{x, Q'^) by Donnachie and Landshoff as a function of for 
the following values of x (from the top): 0.008, 0.035, 0.07, 0.18 and 0.45 together 
with data from NMC 0, SLAC H and BCDMS [||. The F2 values are scaled 
by factors (from the top): 16, 8, 4, 1.5 and 1 for clarity. The error bars represent 
statistical and systematic errors added in quadrature (from 0). 

5. Parametrisation of -^^'(x, Q"^) by Abramowicz et al. The x values, data and scaling 
factors are as in fig. 4 (from 0). 

6. Parametrisation of F2{x,Q'^) by Badelek and Kwiecihski as a function of for 
the following values of x (from the top): 0.008, 0.0125, 0.0175, 0.035. 0.07 and 0.1 



together with data from NMC [H, SLAC H and BCDMS H. The F2 values are 



scaled by factors (from the top): 32, 26, 8, 4, 2 and 1 for clarity. The error bars 
represent statistical and systematic errors added in quadrature (from p). 

7. Regions of the (x, /i^) plane as used in the parametrisation by Schuler and Sjostrand 
(from PI). 



8. Comparison of the parametrisations of F2{x, Q"^) by Abramowicz et al. (full curves), 
Badelek and Kwiecihski (dashed curves) and Donnachie and Landshoff (dash-dotted 
curves) as functions of for the following values of x (from the top): 0.008, 0.02, 
0.05 and 0.1. The values are scaled by the factors (from the top): 27, 9, 3 and 1 for 
clarity (from 0). 

9. Sketch of the x dependence of the ratio F2 / F2 where various nuclear effects are 
indicated. 

10. Multiple scattering of vector mesons. Lines in the lower part of the diagram denote 
nucleons. 

11. Hand-bag diagram for the virtual Compton scattering on a nucleus and its decom- 
position into a multiple scattering series. Lines in the upper parts of diagrams 
denote quarks (antiquarks) and lines in the lower part denote nucleons. 

12. Kinematic ranges covered by different triggers in the NMC (from [p5|). 

13. NMC (preliminary) results - BP function of x compared with the QCD 
predictions (the curve, see |Q for details) and with the results of SLAC (open 



symbols, ||65[). Figure comes from ref. 
14. The F:^ data from NMC compared to the data from SLAC H and BCDMS |5 



The error bars are the quadratic sums of the statistical and systematic errors, ex- 
cluding the overall normalisation uncertainty. The curve is a result of a 15-parameter 
function fit to all three data sets, given in Table 1 (from [ffl]). 



15. Comparison of the (low Q^) ^ dependence of the NMC |l[, EMC NA28 H] and 
SLAC m results (from W). 
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16. Results of the QCD fit to the data. The sohd hne is the result of the QCD fit 
with higher twist included, cf. eq. (^). The dotted curve shows the contribution of 
F2^. In the fit 90 (280) GeV data were renormalised by 0.993 (1.011). The errors 
are statistical (from 0]). 

17. Preliminary measurements of F2{x, Q"^) by the E665 Collaboration at Fermilab. The 
errors are statistical. Systematic errors are 5-15 % (from |[58[| ). 



18. NMC results on the ratio F2/F2 as the function of x at values averaged in 
each X bin. The solid symbols mark the final data set while the open ones are still 
perliminary data taken with a special "small x trigger". Errors are statistical; the 
band at the bottom indicates the preliminary estimate of systematic errors (from 



11) 



19. E665 results for an/cTp as a function of x at values averaged in each x bin for three 
methods of data analysis. Errors are statistical; the total systematic uncertainty is 
less than 3.5% in all x regions. In this figure the NMC data at =4 GeV^ are also 
shown (from ||58|| ). 



20. NMC results on structure function ratios, F^ jF^ as functions of x together with 
the earlier results of SLAC (from [Q). Errors are statistical. 

21. A dependence of the NMC F^ jF^ data for three x bins (from p6|). 

22. Slopes h from the F^ / F^ = a + hlnA fit to the NMC data (from 



23. E665 results on / as a function of x. Errors are statistical, the systematic ones 
are marked as shaded bands. The 3.4% overall normalisation error has not been 
included (from 



24. dependence of the E665 nuclear data in bins of x (from 
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